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We study numerically the two-dimensional flow past a circular cylinder as a prototyp- 
ical transitional flow, and investigate the influence of a generic slip boundary condition 
on the wake dynamics. We show that slip significantly delays the onset of recirculation 
and shedding in the wake behind the cylinder. As expected, the drag on the cylinder 
decreases with slip, with an increased drag sensitivity for large Reynolds numbers. We 
also show that past the critical shedding Reynolds number, slip decreases the vorticity 
intensity in the wake, as well as the lift forces on the cylinder, but increases the shedding 
frequency. We further provide evidence that the shedding transition can be interpreted as 
a critical accumulation of surface vorticity, similarly to related studies on wake instability 
of axisymmetric bodies. Finally, we propose that our results could be used as a passive 
method to infer the effective friction properties of slipping surfaces. 



1. Introduction 

In most continuum studies in fluid mechanics, the goal is to understand and predict the 
behaviour of systems with different physical phenomena at play, but there is in general 
very little debate regarding what the appropriate flow equations and boundary conditions 
are. For gases or liquids flowing under normal conditions in mesoscopic systems (millime- 
ters or larger) bounded by solid walls, it is now universally agreed that the Navier-Stokes 
equations associated with the no-slip boundary condition provide an excellent description 
of the velocity and pressure fields in the fluid. 

Much recent experimental work has been devot ed to the behaviour of fluids at sub- 
mesos copic length scales, microns and nanometers ( Stone et al. ( 2004 ); Squires fc Quake! 



( 20051 )). In the limit where the continuum description remains suitable, the Navier- 



Stokes equations still provide an appropriate predictive framework, but the validity of 
the no-slip bound ary condition has be e n proven experimentally to break down i n ver y 
confined systems (IVinogradoval (|l999h . iGranick. Zhu fc Led (|2003l ). iNeto et al\ (|2005h . 



Lauga. Brenner fc Stond (|2007l )). Wall slip is typically quantified by a slip length, A, 



which is the fictitious distance below the slipping surface where the velocity extrapolates 
to zero. 

Typically, surfaces in contact with a wetting fluid show very little slip, but in a non- 
wetting (or hydrophobic) case, the 'intrinsic' slip length can reach tens of nanometers, 
thereby leading to significant reduction in friction on small length scales. Larger slip 
lengths can be obtained, for example, using super-hydrophobic surfaces, a situation where 
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a high-surface energy liquid in contact with a hydrophobic surface with significant surface 
roughness spontaneously de-wets and transitions from a state where it is everywhere in 
contact with the solid, to a state where it is mostly in con tact wi t h air, and is in contact 
with the solid only at the edge of the surface roughness ( Querel ( 2005f) ). Measurements 
of laminar flow over such surfaces with roughness features in the tens of microns range 
show 'effective' slip lengths of a few tens of mic rons, accompanied by no t ewort hy laminar 
drag reduction (up to 40%, and typically 20%. IOu. Perot fc Rothsteinl (|2004f )). 

Whether the boundary slip is intrinsic to the solid/surface co mbination, or i s use d 
as a generic b u t effe ct ive descript i on of c omplex surface pro cesses ( Vinogradoval ( 19991 ); 
Granick et all (|2003l ). iNeto et all (|2005l ). lLauga et all (|2007l )). the drag reduction it im- 



plies may have a remarkable impact on flows beyond the laminar regime, and in particular 
on the transition to turbulence. This has motivated a few groups to study the effect of 
wall slip on flow at higher Reynolds number. A linear stability study for channel flows 
with slip at the wall revealed that slip delays the occ urrence of flow instabil ity, but has 
virtually no impact on the transient energy growth (jLauga fc Cossul (l2005l) 1. A subse- 
quent numerical study looked at fully-develop ed turbulent chann el flow, analyzing the 
influence of wall slip on turbulent structures ( Min fc Kind |2004)). Surprisingly, while 
slip in the streamwise direction does lead to a decrease in friction, slip in the spanwise 
direction was found to increase turbulent drag. A related work reported computations of 
vortex shedding past a cylinder, both in the laminar and turbulent regimes, with either 
uniform slip alo ng the cylind e r or m ixed slip/no-slip domains distributed along the span 
of the cylinder (|You fc Moinl (|2007l ^. In the laminar regime, a weak enhancement of the 
wake stability was observed, whereas in the turbulent regime a significant decrease of 
the fluctuating forces on the cylinder was noticed, as well as a narrowing of the wake. 
Finally, experiments on rolling droplets wit h Reynolds number in the transitional regime 
showed up to 15% laminar drag reduction (jGogte et all (|2005h ). 

In this paper, we consider two-dimensional flow past a circular cylinder as a prototyp- 
ical transitional flow, and ask the following questions. Generically, how does slip on the 
surface of the cylinder modify the transition from steady to unsteady flow? How does 
it change the stability diagram and the overall flow dynamics? The paper is structured 
as follows. We describe our numerical approach and validate our code in section 2. The 
results of our study, namely the influence of slip on the dynamics of the two-dimensional 
wake behind a cylinder, are discussed in section 3. A possible application is discussed in 
section 4 where we propose to use our results as a passive method to infer the friction 
properties of slipping surfaces. We conclude with some suggestions for future work. 



2. Setup, numerical method and validation 

We consider a two-dimensional incompressible flow past a uniformly slipping circular 
cylinder of radius a. The upstream flow is uniform, of magnitude U . We solve numerically 
the incompressible Navier-Stokes equations for the velocity, u, and pressure, p, namely 

pf^ + u-Vu) =-Vp + 2 j uV-S, V-u = 0, (2.1) 

where p is the fluid density, fi its dynamic viscosity and S = (Vu + (Vu) T )/2 denotes 
the rate-of-strain tensor. Associated with (|2.ip are the impermeability and slip boundary 
conditions on the surface of the cylinder, which read 



n • u = 0, nxu = 2An x (S ■ n), 



(2.2) 



Influence of slip on the dynamics of two-dimensional wakes 



3 



Grid | Parameters | C'd Cl St 



GRID 1 iV r x N e = 80 x 160, 5 = 0.004a 1.348 0.334 0.176 

GRID 2 iV r x N e = 80 x 80, 5 = 0.004a 1.332 0.315 0.166 

GRID 3 N r x N e = 100 x 160, S = 0.002a 1.347 0.331 0.169 

GRID 4 N r x N e = 160 x 160, S = 0.001a 1.339 0.324 0.167 



Table 1. Computational results obtained with three different grids in the case Re = 100, Kn = 
(no-slip cylinder). The drag (resp. lift) coefficients are defined by dividing the time-averaged drag 
force (resp. the maximum of the fluctuating lift force) per unit length by irpa 2 U 2 /2. The Strouhal 
number is defined as St — 2a/ (TU) where T is the vortex shedding period. 



where A is the cylinder slip length and n is the unit normal to the cylinder surface. 
The two dimensionless parameters characterizing the flow are the Reynolds number, 
Re — IpUaj ' ji, and the Knudsen number, Kn = A/a. 

The computations in this paper were carri ed out with the finite volume JA DIM 
code described in deta i l in pr evious studies, e.g. iMagnaudet. Rivero fe Fabre (1995) and 
Calmet fe Magnaudetl (1997). We use a polar grid extending up to r M = 80a. The num- 



ber of nodes is N r = 80 and Ng = 160 along the radial and polar directions respectively. 
A uniform distribution in 9 and a geometrical progression along r are used. To properly 
capture the vorticity generated at the cylinder surface, the thickness 5 of the row of cells 
closest to the cylinder surface is chosen as S = 0.004a, so that at least 7 cells are located 
within the boundary layer for the highest Reynolds number considered, the thickness Sbl 
of the boundary layer being estimated as Sbl ~ aRe~ x / 2 . The influence of these numeri- 
cal parameters was carefully checked to make sure that the results are grid- independent. 
The example reported in Table Q] shows the drag and lift coefficients and the Strouhal 
number obtained for Re = 100 in the no-slip case with the grid used in this study, called 
GRID1. These results are compared with those given by three other grids (see also Table 
(2). It may be observed that the changes in N r , Ng as well as in 5 do not induce significant 
modifications of the drag, lift and shedding frequency (most of the variation on the last 
two quantities are due to the marginal sampling time rather than to the grid). Similar 
results were obtained for the surface vorticity near the rear stagnation point. This region 
is very sensitive to the grid characteristics because the sign of the vorticity changes due 
to separation, and a very good agreement was also observed between the predictions pro- 
vided by the four grids. The cod e has been extensively va lidated in the past, especially 
for bubbles and solid spheres (e.g. IMagnaudet et all (jl995l )1. In order to provide an extra 
validation in the present configuration, the values of the drag coefficient, the maximum 
value of the lift coefficient and the Strouhal number are reported for Kn = (no-slip) in 
Tableland compared with a selection of previous results available in the literature. 

Let us remind that the flow past a no-slip circular cylinder first separat e s for a critical 
Reynolds number Re\ slightly less than 7 ( Chen. Pritchard fc Tavener ( 1995T) ). while 
a second transition where the flow becomes time-dependent and vortices of alternating 
signs are shed from either side of the c ylinder occurs at a nother critical Reynolds number 
Re — Re2 which is between 46 and 50 ((Williamson ( 19961 )1. the Strouhal number St{Re2) 
being about 0.13. Present computations indicate Re x = 6.5, Re 2 = 47.5 and St(Re 2 ) ~ 
0.131, in excellent agreement with previous studies. 

Numerically, the slip condition ()2.2|1 is implemented in the following manner. In a 
system of polar coordinates, the tangential strain rate at the surface is S r g\ r=a = 
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Table 2. Drag coefficient, lift coefficient and Strouhal number as a function of the Reynolds 
number in the case of a no-slip cylinder (Kn = 0) and a shear- free cylinder (Kn = oo). Results 
for the n o-slip cylinder are c ompa re d with previous re s ults fr om D NS (DNS1, DNS2, and DNS3 
are fr om iDennis fc Chang] ||1970D . iPersillon fc Brazal (119981 ) and iKarniadakis fc Triantafvlloul 
l|1989f ). respectively) and experiments l|Williamsonl (|1988T )). 



(dug/dr — u$ I a)\ . Evaluating the Taylor expansion of the tangential velocity in the 
vicinity of r — a at the centre of the first two rows of cells surrounding the cylinder (i.e. at 
distances d\ and c?2 {d\ < (fe) from the surface), and equating these expansions with the 
numerical values found at the previous iteration for u$x = u e\ r=a+ d 1 an d u 92 = u o\ r=a+ d 2 
respectively, yields a second-order accurate explicit expression for dug/dr\ r=a . Combin- 
ing the resulting approximation with the second equation in f|2 . 2|) . the surface velocity 
u o\ r=a i s obtained as 

\di(d 2 -d 1 ) d 2 (d 2 -di) J \X a d x d 2 J 

Using (|2.3p . the surface shear stress ^^S r e\ r=a involved in the local momentum balance 
is finally obtained through the above approximation for due / dr\ r=a . 

3. Influence of slip on the wake dynamics 

3.1. Stability diagram 

The main result of our investigation is illustrated in figure [TJ where we display the 
stability diagram for the wake in the (Re, Kn) plane. The regions separated by the two 
solid lines are the subdomains where the wake is, in the order of increasing Reynolds 
number, steady and unseparated, steady and separated, and unsteady (i.e. in the vortex 
shedding regime). The main conclusion inferred from figure [1] is that slip delays both the 
onset of separation and the onset of vortex shedding. The effect is significant, and small 
values of Kn lead to large changes in the critical Reynolds number Re 2 , with a very 
sharp increase for Kn > 0.5. Note also that the critical line for the onset of separation 
behind the cylinder is not monotonic. That is, for Knudsen numbers typically in the 
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Figure 1. Stability diagram in the (Re,Kn) space and structure of the wake vorticity for 
some selected values of the Reynolds number, Re, and Knudsen number, Kn (color scheme 
is the same for all figures; red is positive vorticity). For each Kn (resp. Re), the transition 
curves were obtained by varying Re (resp. Kn) until the results for two consecutive values were 
found to lie on each side of the corresponding transition; a linear interpolation was then used 
to determine more accurately the corresponding critical Re or Kn. The closed circles on the 
transition curves indicate the pairs (Re, Kn) at which the flow was explored for nearly-critical 
conditions. Separation and vortex shedding were detected by examining the occurrence of a 
change of sign in the surface velocity and that of a nonzero lift force, respectively. 



range [2., 2.5], the flow separates only within a finite range of Reynolds number, say 
[Ren, Rei u ], whereas it remains unseparated both for Re < Ren and for Re > Re\ u (see 
the discussion in 3.3). We also illustrate in figure [IJ as insets, the flow in the wake by 
displaying vorticity snapshots for some selected values of Re and Kn. For a given Re, 
increasing Kn leads to a decrease in the vorticity intensity in the boundary layer and 
wake, and a narrowing of the wake due to the decrease of the separation angle. In the 
case Re — 200, the separation angle varies from 72.6° for Kn = to 27.2° for Kn = 0.5 
and finally falls to zero for Kn = 2.45. 

3.2. Global flow characteristics 

The changes in the forces on the cylinder as a function of wall slip are shown in figure 
fS^-b. In figured, we display the normalized time-averaged drag coefficient C^(Kn) = 
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Figure 2. Drag and lift coefficients as a function of the Knudsen number for Re = 20 
(o, only in (a)), 50 (□), 100 (o), 200 (A), 500 (V) and 800 (*). (a): normalized drag co- 
efficient, Cp(Kn) — (CoiKn) — Cd(oo))/(Cd(0) — Cd(oo)); (b): normalized lift coefficient, 
C* L (Kri) = Cl(Ku) /Cl(0). In (b), the Knudsen number is normalized by its critical value, 
Kn c (Re), above which the flow is no longer unsteady (Kn c =0.09, 0.38, 0.47, 0.49 and 0.53 for 
Re=50, 100, 200, 500 and 800, respectively). 



(Ci)(Kn) — CD(Kn = oo)) /(Cd(Kii = 0) — Cd(oo)). The values corresponding to the no- 
slip limit (Kn = 0) and to perfect slip (Kn = oo) are available in Table 2. As expected, 
the drag on the cylinder decreases monotonically with increasing wall slip, with a steeper 
decrease as the Reynolds number goes up. Typically, the normalized drag coefficient 
crosses the median value 0.5 for Kn « 0.5 when Re — 20 and for Kn ~ 0.1 for Re = 800. 
In other words, the wake becomes more sensitive to a small amount of slip as inertial 
effects increase. We also see on Table 2 that the change in the drag coefficient between the 
no-slip and infinite slip cases is small for small Reynolds numbers, but becomes dramatic 
as the Reynolds number increases. For Reynolds numbers beyond Re = 50, the amount 
of drag reduction is larger than 50%, and is above 80% for Re > 200. In figure [2]d, we 
display the normalized lift coefficient C1(Kn) = Ch(Kn)jCh(Kn = 0), as a function of 
the Knudsen number normalized by the critical value Kn c beyond which the wake no 
longer sheds vortices (and therefore displays no lift). Again, we observe a smooth decay 
of the lift force with increasing wall slip, and an increase of the sensitivity of the wake 
as the Reynolds number is increased. The origin of the drag reduction with increasing 
slip is twofold. First, for a sufficiently large Reynolds number, the vorticity at the surface 
decreases from 0(Re 1 / 2 )U / a for Kn — to 0(l)U/a for Kn = oo (see figure[3}, and the 
skin friction decreases accordingly. Second, this decrease in the surface vorticity reduces 
or even suppresses separation at the cylinder surface. Therefore the fore-aft asymmetry 
of the surface pression distribution decreases, thereby reducing the form drag. As less 
vorticity is produced at the cylinder surface when Kn increases, the strength of the 
vortices shedded in the wake decreases with increasing slip, and so does the lift force. 

Further insight on the wake dynamics can be gained by computing the shedding dom- 
inant frequency 1/T as a function of slip. The results are illustrated in figure [3] where 
we display the normalized Strouhal number for the wake, St*(Kn) = St(Kn)/ St(0), as 
a function of Kn/Kn c for different values of Re. The shedding frequency increases with 
Kn, with up to a factor of 4 increase at large Reynolds number. Hence, although slip 
leads to a decrease of the force amplitudes on the cylinder, it increases the frequency at 
which these force components oscillate. This is because the time rate-of-change of the 
vorticity at the cylinder surface (which is proportional to St) results from both the diffu- 
sion of vorticity across the boundary layer and its advection along the cylinder surface. 
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Figure 3. Normalized Strouhal number St*(Kn) — St(Kn) / St(0), as a function of Kn for 
Re — 50, 100, 200, 500 and 800. St was computed by considering the most energetic peak in the 
lift force spectrum. The definition of symbols and that of Kn c are similar to those in figure [5] 



While the former is almost unaffected by slip, the strength of the latter increases linearly 
with the tangential velocity at the cylinder surface. Combining (|2.2|) with an estimate of 
the boundary layer thickness 5bl/o- = aRe~ x / 2 with a = O(l), one can show that the 
magnitude of the tangential velocity scales as aRe 1 / 2 /(aRe 1 / 2 + 1 + Kn~ x ), so that St 
increases with Kn when Re is kept fixed. 

3.3. Surface vorticity 

There are strong indications that the formation of a standing eddy and the occurrence of 
wake instability past axisymmetric bodies arises when the vorticity accumulated around 
the body, i.e. the difference between the amount of vorticity generated at the body 
surface and that which is evacuated downstream in the wake, exceeds some critical value 
(|Leall (11989^ . This idea explains why the wake past rising spherical gas bubbles remains 



axisymmetric for all Reynolds numbers, while the axial symmetry breaks down when 
the bubble deformation exceeds some critical value, vorticity at a perfect- s lip su rface 
being directly proportional to the surface curvature ( Magnaudet &: Mougin ( 20071 )). In 



the latter paper it was shown that, once expressed in terms of the maximum vorticity at 
the body surface, u) maxi the criterion predicting whether the wake past a solid sphere is 
stable or not is identical to that predicting wake instability for a shear-free bubble. This 
criterion indicates that the wake is unstable when uj max (Re) / (U / a) > g(Re), where g is 
a function that weakly depends on the Reynolds number. In other terms, the strength of 
the vorticity at the body surface is the key quantity that determines the stability of the 
wake, irrespective of the specific dynamic boundary condition at the surface. 

It is of interest to examine how the same argument works in the present configuration. 
For this purpose, we show in figure [4] the quantity uj max a/U, for the two transitions 
displayed in figure [IJ as well as along \so-Kn curves, as a function of the Reynolds 
number. Note that the curve corresponding to Kn = exhibits a 1/2 slope for large 
enough Reynolds number, reflecting the fact that the normalized vorticity at a no-slip 
surface grows like Re 1 / 2 . The results shown in figure Q] reveal that along the shedding 
transition curve, the Reynolds number increases by a factor of 20, whereas the crit- 
ical vorticity, although dependent on i?e, displays little variation and remains in the 
range 5 < Lu max a/U < 7. Hence, present results offer additional evidence that the shed- 
ding transition occurs when the surface vorticity exceeds some critical value that weakly 
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Re 

Figure 4. Maximum dimensionless vorticity Umax /(U/ a) at the cylinder surface as a function 
of the Reynolds number (log-log scale). The thick solid and dash-dotted lines correspond to 
Kn — and Kn — oo, respectively; the thin solid lines correspond to Kn — 0.05, 0.1, 0.2, 0.4, 
1, 2.2 and 5, respectively. The long and short dashed lines delimit the regions where the wake is 
unseparated, steady and separated, and unsteady, respectively. These lines were obtained using 
the procedure described in the caption of figure 1. 



varies with the Reynolds number, irrespective of the dynamic boundary condition at the 
cylinder surface. In particular one can infer from figure 2] that the wake past a shear-free 
cylinder (for which uj max a/U = 4, as shown by assuming a potential flow solution with a 
vortex sheet at the cylinder surface) remains steady and unseparated for arbitrarily large 
Reynolds numbers. The same argument allows us to understand why according to figure 
[TJ the wake separates only within a finite range of Reynolds number, [Ren, Rei u ], when 
Kn exceeds some critical value: vorticity generation at the surface increases only weakly 
with Re in presence of significant slip, while the vorticity flux advected downstream in 
the wake increases linearly with the upstream velocity (and therefore with Re), imply- 
ing the existence of an upper value of Re beyond which there is not enough vorticity 
accumulated behind the cylinder for the wake to separate ( Leal ( 1989h ). 



4. Application to the effective friction of slipping surfaces 

The results presented in figureQ]could be used as a method to passively infer the friction 
properties of surfaces prone to slip. For example, given the potentially large drag reduc- 
tion resulting from super-hydrophobic coatings, there is growing interest in obtaining 
quantitative information regarding their overall friction properties, and the effectiveness 
of particular geometrical designs (length scales and morphology). The methods used so 
far in the literature involve real friction measurements, where a particular surface is mod- 
ified to become super- hydrophobic, a p ressure gradient is applied, and the change in the 
flow rate is measured ( Ou et all ( 20041 )). These are arguably difficult experiments and it 
would be useful that simpler methods be available. 

In this spirit, we propose to coat a cylinder with the studied slipping surface, and to 
place it in a uniform flow field with known upstream velocity U . As the flow speed is 
increased beyond a critical value, the flow will become unsteady, and we can then exploit 
our numerical results to infer the effective slip length on the surface, without performing 
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any actual friction measurement. With the knowledge of the critical value of U, and 
hence Re, at which vortex shedding sets in, one reads the corresponding value of Kn 
from figure [TJ and infers the effective slip length describing the surface. According to the 
results displayed in figure [TJ this will work as long as Kn < 0.5. 

For super-hydrophobic surfaces, this method will be effective provided two hypotheses 
are satisfied. We first need to ensure that the surface remains super-hydrophobic, and 
does not transition back to a wetted state. We further need to ensure that the effect 
be sufficiently strong to be experimentally measurable. Regarding the first point, it is 
necessary that the flow-induced pressures, which are on the order of pU 2 , are small enough 
to prevent transition to the wetting state. For typical super-hydrophobic surfaces made of 
a dense population of pillars with size r, the critical pressure is a ca pillary pressur e , on t he 
order of o~/r, where a is the surface tension of the liquid in air ((Bartolo et all (|2006h ). 



Wetting will therefore be avoided if pU 2 < a/r, or Re < a/y/Ir, where £ — n 2 /pcr 
is the Ohnesorge length for the liquid. Regarding the second point, we need the slip 
len gth to be suf f icientl y large for the effect to be detectable. According to previous work 
bv lYbert efaH (|2007t ). we expect A ~ d, where d is the distance between the pillars 



composing the super-hydrophobic surface. To fix ideas, if we wish the transition to occur 
at a particular Kn, we need a as d/ Kn. This limit will be compatible with the constraint 
found above if ReKn < d/y/Ir and, with r ~ d, we need therefore to be in the limit 
where (Re Kn) 2 < r/£. For an air- water system, we have I as 14 nm. The critical Reynolds 
number for, say, Kn — 0.1 is about 60, which is sufficiently above Re2(Kn = 0) = 47.5 
to be measured without ambiguity, so r and d need to be at least a few microns, and the 
proposed experiment can be achieved in a microfluidics device. 



5. Concluding remarks 

In this work, we performed direct numerical simulations to study the influence of 
generic slip boundary conditions on the dynamics of the two-dimensional wake behind a 
circular cylinder. Specifically, we have shown that slip: (1) delays the onset of recirculation 
and vortex shedding in the wake; (2) decreases both drag and lift; (3) decreases the 
vorticity intensity in the near wake; (4) increases the shedding frequency. We have also 
discussed a practical application of our results and proposed that they could be used as 
a passive method to infer the effective friction properties of slipping surfaces. 

The most severe limitation of the current work is the two-dimensional assumption 
for all Reynolds numbers. In the case of a long, no-slip cylinder, three-dimensional ef- 
fects are known t o become important when the Reynolds number reaches Re as 190 
IWilliamsonl(|l996t ). However, since these three-dimensional effects occur well into the un- 
steady regime, the wake is likely to remain two-dimensional close to the shedding tran- 
sition, and our critical curves Re\ = f(Kn) and Re-i = g(Kn) should not be affected by 
these effects. Regarding the applicability of our results as a method to passively estimate 
wall slip, other effects such as geometrical confinement would lead to changes in the de- 
tails of the results shown in figure [TJ and each experimental setup would be characterized 
by its own transition curve to be obtained numerically, but the general principle remains 
valid. Moreover, in general, most complex surfaces cannot be adequately modeled by a 
single homogeneous slip coefficients, but instead a detailed look at the geometrical and 
physico-chemical nature of the interface is required. Results in the case of micro-textured 
surface will be reported in future work. 
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